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The effects of nonmagnetic impurities on tfie Q = (tt, tt) spin-fluctuation spectral weight fmx(Q, i^) 
are studied within the framework of the two-dimensional Hubbard model using the random phase 
approximation. In the first part of the paper, the effects of the nonmagnetic impurities on the 
magnetic susceptibility of the noninteracting {U — 0) system, xo(q, <^), are calculated with the self- 
energy and the vertex corrections using various forms of the effective electron-impurity interaction. 
Here, the range and the strength of the impurity potential are varied as well as the concentration 
of the impurities. It is shown that the main effect of dilute impurities on xo(Q,'jj) is to cause a 
weak smearing. In the second part, Imx(Q,<^) is calculated for the interacting system. Here, the 
calculations are concentrated on the processes which involve the impurity scattering of the spin 
fluctuations with finite momentum transfers. In order to make comparisons with the experimental 
data on the frequency dependence of Imx(Q,w) in Zn substituted YBa2Cu307, results are given 
for various values of the model parameters. 

PACS Numbers: 74.62. Dh, 74.72. Bk, 75.10.Lp, 74.25.-q 



I. INTRODUCTION AND THE MODEL 

Zn substitution has been used as a probe of the elec- 
tronic properties in the layered cuprates and has yielded 
important results. Small amount of Zn impurities lead 
to a strong suppression of d-wave pairing Resistiv- 
ity measurements find that Zn impurities act as strong 
scatterers Local magnetic measurements on Zn sub- 
stituted cuprates have also yielded valuable information 
It has been found that within the presence of Zn 
impurities the uniform susceptibility of YBa2Cu307_5 
follows a Curie-like temperature dependence. 

Inelastic neutron scattering experiments find that Zn 
impurities cause important changes in the Q = (tt, tt) 
spin-fluctuation spectrum ^-Q. The experiments have 
been carried out for 2% [|| and 0.5% j|] Zn concentra- 
tions. It is found that even 0.5% Zn substitution in 
YBa2Cu307 induces a peak in lmx{Q,i^) at w = 40 
meV in the normal state The uj width of the peak 
is about 10 meV and its width in momentum space is 
resolution limited. This peak becomes observable below 
250° K and its intensity increases as T is lowered with- 
out a discontinuity in its slope at Tc — 87°K. These re- 
sults are quite different than those on pure YBa2Cu307, 
where the peak is observed only in the superconducting 
state |8-10 . While in the case of 0.5% Zn substitution, 
no spectral weight is observed below ~ 35 meV, in the 
2% Zn substituted sample there is significant amount of 
spectral weight down to 5 meV. In this case, there is also 
a broad peak at 35 meV. The momentum width of this 
peak is 0.5A^^, which is about twice that of the peak in 
the 0.5% Zn substituted and the pure samples. Hence, 
the two main features induced in Imx(Q,w) in the nor- 
mal state are the peak and the low frequency spectral 
weight observed in the 2% Zn substituted sample. Both 



of these features become observable already above Tc- 
The calculations reported in this article were carried out 
in order to investigate the origin of these two effects on 
the frequency dependence of Imx(Q,u;) in YBa2Cu307 
in the normal state. 

Here, the results of diagrammatic calculations on the 
effects of nonmagnetic impurities on Imx(Q, i^) obtained 
using the framework of the two-dimensional Hubbard 
model will be given. In the first part of the paper, the 
effects of various types of effective electron-impurity in- 
teractions on the magnetic susceptibility of the nonin- 
teracting {U — 0) system, xo(Q7'j-')i will be calculated. 
In this case, both the strength and the range of the 
impurity interaction, in addition to the concentration 
of impurities, will be varied. It will be found that in 
all these cases the main effect of impurity scattering is 
to cause a weak smearing of the structure in xoCQi'^)- 
Hence, at this level, if xo is used in an RPA expression, 
x(Q,w) = Xo(Q,w)/(l-f/xo(Q,w)), then one obtains a 
smearing of Im x(Q, by the impurity scattering rather 
than an enhancement as seen in the experiments. Here, 
it will be also noted that, for 2% impurity concentration, 
scatterings from an extended impurity potential could 
induce spectral weight at low frequencies. 

In the next part of the paper, the effects of the pro- 
cesses where the spin fluctuations are scattered by the im- 
purities with finite momentum transfers are calculated. 
This type of processes will be called the "umklapp" pro- 
cesses, since they involve the scattering of the spin fluc- 
tuations by the impurity potential with finite momen- 
tum transfers as in the case of the scattering of the spin 
fluctuations by a charge-density-wave field. In this part, 
the effects of the umklapp scatterings will be estimated 
by calculating the irreducible off-diagonal susceptibility 
Xo(q, q', ijj), where q q', in the lowest order in the 
strength of the impurity potential. It will be seen that the 



1 



important umklapp processes are the ones which involve 
the transfer of momentum 2k.p to the spin fluctuations, 
and that they lead to a peak in Imx(Q, i^) at ujq — 2|^|, 
where fx is the chemical potential. The underlying rea- 
son for this is a kinematic constraint which prohibits the 
creation of a particle-hole pair with center of mass mo- 
mentum q = (tt, tt) — 2\ip and energy lo > ujq. This con- 
straint causes a nearly singular structure in xoiQ^I,^) 
at uj ~ luq, which in turn leads to the peak in Imx(Q, i^) 
at uq. This effect has been noted previously Here, 
results will be given for various sets of the model parame- 
ters in order to make comparisons with the experimental 
data. While in the pure system at low temperatures, 
the Q = (7r,7r) spin fluctuations are gapped below luq, 
through this type of scatterings the Q = (tt, tt) spin fluc- 
tuations get mixed with the other wave vector compo- 
nents which are not gapped. This process can also lead to 
finite spectral weight for lu < ujq. Comparisons with the 
experimental data suggest that the umklapp processes 
play an important role in determining the Q = (tt, tt) 
spin dynamics in Zn substituted YBa2Cu307. 

The starting point is the two-dimensional single-band 
Hubbard model given by 



-f^^4a(^^-y^ (1) 



which will be used to model the spin fluctuations of the 
pure system. Here Cia (cj^.) annihilates (creates) an elec- 
tron with spin a at site i, t is the near-neighbor hopping 
matrix element and U is the onsite Coulomb repulsion. 
For simplicity, the hopping i, the lattice constant a and 
h are set to 1. 

Within the presence of an impurity, the magnetic sus- 
ceptibility is defined by 



X(q,q', jw„ 



dre" 



'(m-(q,r)TO+(q',0)), (2) 



where m+(q) = N Ep Cp+q|Cpj., m (q) = 
(m"'"(q))^, and ujm ~ 2m'KT. By letting iuim ^ uj + iS, 
one obtains x(q, q',(jj). If one assumes that the effec- 
tive interaction between an impurity and the electrons 
can be approximated by a static potential, then the RPA 
expression for the magnetic susceptibility becomes 



x(q, q',^)= xo(q, q' 



xo(q, q", ui)xic^', q', t^), 



(3) 



where xo(qi q': uj) is the irreducible susceptibility dressed 
with the impurity scatterings. The off-diagonal q 7^ q' 
terms of xo{^., q'l vanish for the pure system, but they 
are finite within the presence of impurity scattering. Af- 
ter obtaining x(q, q', ui) from Eq. (0), the averaging over 



the impurity location can be done, which sets q = q'. 
x(q, q, w) obtained this way neglects the interactions be- 
tween the impurities, and here it will be assumed that 
these can be ignored in the dilute limit. In the following, 
for simplicity, x(q, q, and xoil^ljUj) will be denoted 
by x(q,w) and xo(q,w), respectively. 

In solving Eq. (Q), if the off-diagonal terms are omitted, 
then one obtains 



x(q,^^) 



xo(q:^) 
1 - c/xo(q,' 



(4) 



where xo(qi^) is dressed with the impurity scatterings. 
In the next section, xo(Q = {'^,t^),uj) will be calcu- 
lated using various impurity potentials. In the third sec- 
tion, x(Q,'jj) will be calculated without omitting the off- 
diagonal terms in Eq. (0). 



II. EFFECTS OF THE NONMAGNETIC 
IMPURITIES WITHOUT THE UMKLAPP 
SCATTERINGS 

In this part, the effects of dilute nonmagnetic impuri- 
ties on the magnetic susceptibility of the noninteracting 
system will be calculated using various impurity poten- 
tials. The method used here for calculating xo(Q,w) is 
similar to those used in Refs. [p"2|-p^. Both the self- 
energy and the vertex corrections induced by the impu- 
rity scattering will be included | p^ , p^ . 

In Refs. [ p7[p^ ], it has been shown how the electronic 
correlations lead to an extended effective interaction be- 
tween an impurity and the electrons. The extended na- 
ture of the effective impurity potential has been also em- 
phasized in Ref. p9| . Assuming that it can be approxi- 
mated by a static form, the potential due to an impurity 
at site Tq can be written as 



eff 



(5) 



where u denotes the distance from the impurity and a de- 
notes the different partial wave components. The single- 
particle operators c^aa are given by 



(6) 



where sums over the sites at a distance v away from 
the impurity and g^aiduYs are the coefficients of the 
partial- wave components. In the following, an impurity 
interaction with a range of lattice spacings, which in- 
cludes the second near-neighbor site, will be considered. 
For simplicity, subscript i will be used to denote both v 
and a. 

Within the presence of impurity scattering, the single- 
particle Green's function defined by 
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G{v,r) = -{TrC^.{T)clM) (7) 
is obtained from 

G{p,iujn)= Go(p,it^„) 

+n,(Go(p,icj„))2^5,(p)Ty(it^„)5j(p), (8) 

ij 

where is the concentration of the impurities and the 
Matsubara frequency LOn = (2n + l)7rT. For an impu- 
rity potential with a range of \/2 lattice spacings, the 
subscripts i and j vary from 1 to 9. The single-particle 
Green's function of the pure system Go(p, jti-'n) entering 
Eq. (0) is given by 



Go(p,«a;„) 



1 



ILOn - e-r 



(9) 



where the single-particle dispersion relation is 

Sp — —2t{cos px + cos py) — ^. (10) 

The terms contributing to G(p,ia;„) are illustrated di- 
agrammatically in Fig. 1(a). In Eq. (P, Tij{iujn) is the 
impurity-scattering i- matrix and gi (p) 's are the form fac- 
tors which are given by 



(11) 



where i denotes {v,Oi). The first five of the nine form 
factors used here are 

ffOs(p)= 1 

gis(p)= cospx + cospy 
ffipx(p)= sin 
5ip„(p)= sinpy 

5id(p)= cospx - cospy (12) 

with similar expressions for the remaining a = com- 
ponents having s, px, Py and d-wave symmetries. The 
t-matrix Tij{iLUn) is obtained by solving 

Tij{iujn) = S^jVi + ^ Fii{iuj„)Tij{iujn), (13) 



where 



F^j{iujn) = -^^gj(p)Go(p,it^„)5j(p). 



(14) 



In order to calculate xo(Qj*i^m), the irreducible inter- 
action in the particle-hole channel due to impurity 
scattering is needed. In Fig. 1(b), F'' is illustrated dia- 
grammatically and the corresponding expression is 

r°'jy(P,p',*^«,«w„) = -nig,(p)5i-(p + Q) 

xT^j{iujn)Ti'j>{iUn + iujm) 9 ]{p') 9 rip' + Q)- (15) 
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FIG. 1. (a) Feynman diagrams for the single-particle 
Green's function dressed with the impurity scatterings, (b) 
Irreducible vertex P" in the particle-hole channel arising from 
the impurity scattering, (c) Reducible vertex F in the parti- 
cle-hole channel given in terms of F". (d) Irreducible magnetic 
susceptibility xo(Q, j<^m) dressed with the impurity-induced 
self-energy and vertex corrections. 

Here, F*' is evaluated with the center-of-mass momentum 
Q = (tt, tt). Next, the reducible vertex in the particle- 
hole channel is obtained by solving 

Fii'jj'(p,p',itj„,ia;m) = T%, jj, {p,p' ,iLUn,iLu„i) 



^ — ^ F° ,_„,(p, k, iuJn,iuj,n)G(k, iuj^) 

U,ss' k 

xG(k + Q,icj„ +ia;„)Fss/j.,v(k, p',zcj„,zwm), (16) 

which is illustrated in Fig. 1(c). Here, it is noted that F° 
is calculated using the bare single-particle Green's func- 
tion Gq, while F is calculated using the Green's function 
dressed with the impurity scatterings. This is necessary 
in order to prevent double-counting p^ ] . In terms of the 
reducible vertex, XoiQ^i^m) is given by 



IUJ„ 



X 2^ Tu>jj>{p,p',iujn,iuJm)G{p' + Q,iuJn+iu}m.) 

xG(p',zc^„). (17) 

Here, Xo(Qi ^'^m) includes only the impurity induced self- 
energy corrections, and it is given by 

Xo(Q, ii^m) = X! + + 



xG(p,iw„). (18) 
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These results in terms of the Matsubara frequencies are 
analytically continued to the real frequency axis by the 
Pads approximation. In the following, the results on 
Xo(Q,^) obtained this way will be compared with the 
Lindhard susceptibility of the pure system, 



v^rn , A - /(gp+q) - /(gp) 



(19) 



Figure 2 shows results obtained using a strongly at- 
tractive onsite impurity potential Vq = —20 and an im- 
purity concentration of rij = 0.02. In addition, here fill- 
ing (n) = 0.875 and temperature T = 0.02 arc used. In 
Fig. 2(a), XoiQji'^m) versus LOm is shown. Also shown 
in this figure are Xo(Qi^'^m), which does not include the 
impurity vertex corrections, and the Lindhard suscepti- 
bility Xo {Q,i^m) of the pure system. Figures 2(b) and 
(c) show the corresponding real and imaginary parts ob- 
tained by the Fade analytic continuation. Here, one ob- 
serves that the impurity-induced self-energy corrections 
suppress Xo(Q) if^m) and smear the structure in xo(Q7 w). 
For instance, the hump in Re Xq(Q, cj) at a; sa 2|/i| ~ 0.48 
is smeared by the self-energy corrections, but when the 
vertex corrections are included, the efi^ect of the 
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FIG. 2. Results on (a) xo(Q,iWm) versus Wm, (b) 

Rexo(Q,'i^) versus and (c) Imxo(Q,<j.') versus to for an on- 
site impurity potential with Vb = —20 and impurity concen- 
tration jii = 0.02. These results were obtained for (n) = 0.875 
and T = 0.02. 



self-energy corrections is nearly canceled. This hump is 
because of a logarithmic singularity in Rexo(Q)<^) at 
T = originating from the dynamic nesting of the Fermi 
surface. In a real system, the deviations from the simple 
tight-binding model which has only near-neighbor hop- 
pings would lead to a suppression of this hump. In ad- 
dition, the scattering of the single-particle excitations by 
the spin-fluctuations would suppress it also. In Figs. 2(b) 
and (c), the difference between the solid and the dot- 
ted lines is of order the resolution of the Fade analytic 
continuation. These calculations were repeated using 
Ui = 0.005 instead of n, = 0.02, in which case the dif- 
ference between xo and Xo becomes even smaller (not 
shown here). 

Figure 3 shows results also for = 0.02 and an onsite 

impurity potential but now with Vq = — 1. One observes 
that in this case the impurities have a stronger effect on 
Xo- In addition, it is noted that Rexo(Q)W) for oj ^ 0.5 
is suppressed, and Rexo(Q,t^ ^ 0) gets enhanced by a 
small amount. Among the various forms tried for the 
impurity potential, this is the only case where an en- 
hancement of xo(Qi'^) by the impurity scatterings has 
been obtained. This effect is due to the enhancement 
of the single-particle density of states at the impurity 
site by the attractive potential. The enhancement of 
Rexo(Q)W ~ 0) at small frequencies could lead to some 
enhancement of the low-frequency antiferromagnetic spin 
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FIG. 3. Results similar to those in Fig. 2 but for Vo = —1. 



4 



fluctuations for a system with a large Stoner factor. How- 
ever, it is seen that Rexo(Q;i^) for ^ near coq = 2\^\ is 
suppressed, and actually it would not be possible to ex- 
plain the peak in the neutron scattering data with these 
results. This calculation was repeated using Vq = 1, and 
a suppresion of Rexo(Q) was obtained due to the de- 
pletion of the single-particle density of states (not shown 
here) . 

Figure 4 shows results on xo obtained using an ex- 
tended impurity potential for rii = 0.02 and 0.005. These 
results were obtained for a potential with a range of 
•\/2 lattice spacings and with the following parameters: 
Vos = -20, Via, = 0.5 and = -0.25 where a de- 

notes the s, Px, Py and d-wave components. These values 
for Vjya's are comparable to those obtained in Ref. 
The calculations were repeated using various other values 
for the Vjya's, and it has been found that small changes in 
Fi/q's do not change the results shown here significantly. 
For instance, increasing V^a's by 50% does not change 
the conclusions of this section. In Fig. 4, it is seen that 
for Tii = 0.02 an extended potential leads to significant 
smearing of the structure in xo{Qt^)- In this case, the 
hump in Hexoi^,'-^) is rounded off, and spectral weight 
is induced for lu < luq- Hence, comparing with Fig. 2, 
one observes that while an onsite impurity potential 
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FIG. 4. Results similar to those in Fig. 2 but for 
an extended potential with the following set of parameters: 
Vos — —20, Via — 0.5 and — —0.5. Here results are 

given for both rii — 0.02 and 0.005. 



does not lead to spectral weight for lu < ujq, an extended 
potential induces spectral weight in the gap. Also shown 
in Fig. 4 are the results for Ui = 0.005, in which case 
the effect of the impurities on xo{Q,^) is weaker. The 
fact that 0.5% impurities induce less spectral weight in 
the gap compared to the 2% case is consistent with the 
neutron scattering data by Refs. ||,^. However, even 
for Ui = 0.02, the amount of the spectral weight in- 
duced in the gap is small. In the next section, it will 
be seen that the umklapp processes could also contribute 
to Imx(Q,aj) for w < loq. 

From the results presented here, one observes that if 
Xo{Q,t^) computed in this section is used in Eq. (||), 
which omits the umklapp scatterings, then one would 
only obtain a smearing of Imx(Q,Ci;) by the impurities. 
Hence, it would not be possible to explain how 0.5% Zn 
impurities induce a peak in Imx(Q,a;) in the normal 
state of YBa2Cu307. In the next section, the effects of 
the impurity scatterings with finite momentum transfers 
will be taken into account. 



III. EFFECTS OF THE IMPURITY INDUCED 
UMKLAPP SCATTERINGS 

Here, Imx(Q,cij) will be calculated without omitting 
the off-diagonal components in Eq. (j^). These off- 
diagonal terms will be calculated in the lowest order in 
the strength of the impurity potential, as illustrated di- 
agrammatically in Fig. 5(a). In the previous section, it 
was found that 0.5% and 2% impurities cause only a weak 
smearing of xoiQ,^)- For this reason, the diagonal com- 
ponents of Eq. (0) will be approximated by the Lindhard 
susceptibility Xq- This will not change the nearly singu- 
lar contribution originating from the umklapp scatterings 
at w luq. The expression for xo(Q; Qj i^m) correspond- 
ing to the diagrams shown in Fig. 5(a) is 



Xo(Q,q,«w™) = -Vo^ X! |G'o(p + Q,«t^n- 

xGo(p,iw„)Go(p + Q - q, iw„) 
-fGo(p, iw„)Go(p + q,«w„ +iujm) 



ILOm) 



xGo(p + Q,iuJn + i^^ni) 



(20) 



Here Q = (tt, tt) and q = Q — Q*, where Q* is the 
momentum transferred during the scattering from the 
impurity. Upon carrying out the summation over iuon in 
Eq. (EG) and letting iujm uj + iS, one obtains 



Xo(Q,q,w)= A(Q,q,w) 

+^(Q,q,^) 



-A*(Q,q,-w) 
B*(Q,q, -uj), 



(21) 



where "*" stands for complex conjugation, and A and B 
are given by 
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A(Q,q,cj) 



^ ^ I ^ - (£p+Q - £p+Q* ) + i& 



(ep+Q - £p) +iS j Sp+Q- - £ 



(22) 



5(Q,q,^) 



TV ^ 



1 



£p+Q' — £p 



/(£p+Q* 



(cj - (£p+Q - £p+Q.) + i6){uj - {sp+Q - ep) + iS) ' 



(23) 



A general form for the effective interaction between the 
electrons and an impurity located at site Tq is given by 



p+kcr'-pc- 



(24) 



p,<T 



Below, it will be seen that for Q* near 2ki?, xo(Q5qi^) 
has a nearly singular structure at w « wq; while for Q* 
away from 2ki?, xolQjQj'j-') is a smooth function of lu 
with a small amplitude. Hence, x(Qi^^) will be calcu- 
lated using only the ±Q* components of the effective 
impurity interaction. 



^0 I](4+Q*'^^P'^ + Cp-Q.^cp^ 

pr- 



(25) 



where Vq = Vq» is taken as a parameter. This is neces- 
sary, since, in order to have sufficient frequency resolu- 
tion, the calculation needs to be carried out on a large lat- 
tice, which is difficult to do using directly Eq. (p4|). Fur- 
thermore, the detailed k dependence of Vk is not known, 
especially for k near 2kF. The scattering of a quasiparti- 
cle with Q* = 2ki? momentum transfer is sketched in the 
Brillouin zone in Fig. 5(b). Using the interaction given 
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FIG. 5. (a) Feynman diagrams for the lowest-order 
terms contributing to the off-diagonal irreducible susceptibil- 
ity Xo(Q, q, ^\ (b) Sketch in the Brillouin zone of a quasipar- 
ticle scattering process where momentum 2k_F is transferred. 



in Eq. (g5|), one obtains for x(Q,w), 
X(Q,C^) = {x^(Q,c^)(f - t/x^(q,c^)) 

+4C/(xo(Q,q,c^))'} X {(I - C/x^(Q,c^)) 

x(l - t/x^(q,c^)) - 4((7xo(Q,q,c^))'}"\ (26) 

where q — Q — Q*. Here, the factor of 4 multiply- 
ing (xo(Q, q, t^))^ is to take into account the scatter- 
ings with momentum transfers (±(5*,T'3*) in addition 
to (±g*,±Q*), where Q* = (Q*,Q*). 

In the following, results will be shown for (n) = 0.86 
and T = 0.05, in which case loq = 2|/i| w 0.55. Fig- 
ure 6(a) shows the real and the imaginary parts of 
Xo(Q,q, w) for Q* = 2ki?. Here, the Fermi wave vector 
ki? has been taken along (1,1) for simplicity. In evaluat- 
ing x(Q,^) with Eq. (|26|), the Lindhard susceptibilities 
Xo(Q,t^) and Xo(q = Q ~ Q*i^) are also used and, 
hence, they are plotted in Figures 6(b) and (c). Here, 
one notes the similarity between the lo dependence of 
Xo(Q,q, w) and Xo(q = Q ^ Q*:'^)- Both have van- 
ishing spectral weight for cj > wg because of kinematic 
constraints. 

In Fig. 7, Imx(Q,'^) versus w obtained from Eq. ( |2^ ) 
by using the results of Fig. 6 are shown for different values 
of U . The solid lines are for Vq = 0.05 and the dashed 
lines are for Vq = corresponding to the pure case. Here, 
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These results were obtained for (n) = 0.86 and T = 0.05. 
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and (c) U = 2.0. These results were for Vo = 0.05 and 0. 
value of Vo is increased to 0.06, the sharp peak at ujo f 
diverges. 



If the 
i 0.55 



it is seen that a peak is induced at coq ~ 0.55 by turning 
on Vo- One also observes that as U is increased, a hump 
develops below the peak. This is because of the RPA 
enhancement of x in the pure case. In a real system, it 
is expected that this hump will be smaller because of the 
band-structure effects and the damping of the quasipar- 
ticles by the spin fluctuations. In Fig. 8, it is also seen 
that the impurity contribution to the low frequency part 
of Imx(Q,a;) increases as U is increased. In addition, 
in Fig. 6(c) it was seen that xo(q = Q — 2kF,tf) has 
sharp structure near wq, but this is not responsible for 
producing the peak in Imx(Q,w). For instance, taking 
out the factor of (1 — Uxoil,^)) from both the numera- 



tor and the denominator in Eq. (|2g) does not change the 
structure of Imx(Q, w). 

Next, in order to have a better understanding of these 
results, a sketch of the important wave vectors and fre- 
quencies are given in the q-w plane in Fig. 8. Here, 
Q = (7r,7r), and q and 2]<.p are taken along (1,1). The 
vertical dashed line denotes q — Q — 2]<.p and the hor- 
izontal dashed line is for uJo ~ The shaded area 
represents the region where Imxo(qi^) 7^ at T = 0. 
Here, one observes that the Q* = 2ki? scattering of the 
Q = (7r,7r) spin fluctuations will lead to a mixing with 
the q = Q — 2ki? component of the spin fluctuations, 
which has spectral weight only for uj < coq. A general 




Q-2kp 



2kp Q 



FIG. 8. Sketch in q — tj plane of the important wave 
vectors and frequencies (not to scale). The vertical dashed 
line denotes q — Q — 2kF where Q = (tt, tt) and the horizontal 
dashed line denotes uo ~ 2|/i|. The shaded area represents 
the region where lmxo{(l,ijj) 7^ at T = 0. Here q and 2k_F 
were taken along (1, 1). 

impurity potential such as Eq. ( p4[ ) would lead to a mix- 
ing with all wave vectors. Hence, the uniklapp processes 
could also contribute to Imx(Q,w) for < ojo, in addi- 
tion to the results seen in the previous section. 

In order to show that the nearly singular uj dependence 
of xo(Q, q, ^) occurs only for Q* sa 2ki?, in Fig. 9 results 
on xo(Q, q, ^)/Vo are shown for various values of Q*. As 
seen in Fig. 9(a), for Q* = 1.02(2ki?), there is a sharp 
structure in xo(Q) Qj ^) ^ less than wqi while for Q* — 
0.98(2ki;-), this occurs at > cjq. As Q* moves away 
from 2ki?, the position of the structure in xo(Qjqi^) 
shifts away from and its amplitude decreases. In 
Figs. 9(b) and (c), results on Rexo(Q: ^)/^ versus 
w are shown for Q* along (1,1) and (1,0), respectively. 
Here, Imxo(Q- ^) /^o is not shown since it is a smooth 
function of lo with amplitude less than 0.05. In these 
figures, it is seen that the magnitude of xo(Q,q, w) at 
w w Wo is considerably smaller when Q* is away from 
2ki?. This supports the use of only the Q* component of 
Vefj in solving for x(Q,'^). If a general potential such 
as Eq. (|2j) instead of Eq. (|2^) were used in solving for 
x(Q,tf), then the peak in Imx(Q,ti^) would again occur 
at Wo but with a broadened width because of the con- 
tributions originating from scatterings with Q* slightly 
away from 2ki?. 

In Fig. 7, one notes that the line shape of the peak 
depends on the value of U . For U = 1.5 and 1.7, the 
peak is asymmetric; there is a dip below the peak. For 
U = 2.0, the dip is not observed. In order to have a 
better understanding of this, further results on the line 
shape are shown in Figures 10(a) and (b) for U ~ 1.7 
and 2.0, respectively. In these figures, if the value of Vq 
is increased to 0.06, the peak at ujq diverges. Also shown 
here in Fig. 10(c) is the quantity 
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FIG. 9. (a) xo(Q, q, 'j.')/Vb versus uj for different values 
of Q* near 2'k.p along (1,1). Here, the numbers next to the 
curves indicate the magnitude of Q* in units of 2fcF • The solid 
and the dotted curves represent the real and the imaginary 
parts. In (b) and (c), results on Re xo(Q, Q, '^)/l^ versus uj 
are plotted for Q* away from 2kF. 



7H=4C/2(xo(Q,q,c^))' 



(27) 



which enters Eq. (p6[). For U — 1.7, the system is away 
from the magnetic instability and, in this case. Re 7(0;) 
has a stronger effect than the imaginary part. While for 
UJ > loq, Re7(w) leads to a stronger RPA enhancement 
of Imx(Q,w), for w < Wo it suppresses luixiQ^uj). On 
the other hand, for [7 = 2, where 1 — URexo{Q,uj) is 
small, the pure system is already close to the magnetic 
instability and Ulmxo{Q,uj) acts as a damping of the 
RPA enhancement. In this case, Im7(ti;) becomes more 
important and it leads to the peak in Imx(Q,w) by sup- 
pressing the damping of the RPA enhancement. Possibly, 
this is the case more applicable to YBa2Cu307. 

The results seen in Fig. 10 were obtained using a finite 
broadening S = 0.001 in Eqs. (H) and (H). For 

comparison, results obtained using 6 = 0.01 are shown 
in Fig. 11. In this case, the structure in 7(0^) and the 
peak in Imx(Q,w) are broader. It is also seen that the 
integrated spectral weight in the peak increases with 6 
for U — 2. If ^ is increased further, the width of the peak 
in Imx(Q, uj) for U = 2 continues to increase (not shown 
here) . It is noted that the scattering of the quasiparticles 
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FIG. 10. Imx(Q,'^) versus uj for different values of Vb, 
and for (a) U = 1.7 and (b) U = 2. (c) ^(uj) versus uj for 
Vb = 0.05 and U = 2.0. These results were obtained for 
S = 0.001. 

by the spin fluctuations could have an effect similar to 
that of the finite broadening S used here. This could be 
another reason for why the peak observed in Imx(Q,w) 
of Zn substituted samples has a finite width. So, in these 
figures it is observed that the quantitative features of the 
changes induced by the umklapp processes depend on the 
model parameters. 



IV. DISCUSSION 

In summary, the effects of dilute nonmagnetic impu- 
rities on the frequency dependence of ImxiQ, uj) have 
been studied. These calculations were motivated by the 
neutron scattering data of Refs. on 2% and 0.5% Zn 
substituted YBa2Cu307. The origin of the peak and of 
the low-frequency spectral weight induced by the impuri- 
ties are the two issues which have been addressed in this 
paper. 

In pure YBa2Cu307, a resonant peak is observed only 
in the superconducting state at w = 41 meV |]8|-|l0|. The 
width of this peak is resolution limited as opposed to 
that observed in the Zn substituted samples. The the- 
ories attribute the resonant peak of the pure sample to 
instabilities in the particle-particle [^U] or the mag- 
netic channels |22-p4| . The calculations presented here 
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FIG. 11. Results similar to those shown in Fig. 10 but for 
(5 = 0.01 and different values of Vo- 

show how a peak in Imx(Q,a;) of Zn substituted 
YBa2Cu307 could arise from the magnetic channel in 
the normal state. However, clearly, the contributions to 
the peak from other channels are not ruled out. 

The calculations have been carried out first without 
taking into account the umklapp scattering of the spin 
fluctuations. Here, the influence of the range of the im- 
purity potential on xo(Q;^) has been studied. When 
a strongly attractive onsite impurity potential is used, 
it has been found that 2% impurities have a negligible 
effect on Xo(Q)'a')- This is in agreement with the cal- 
culations of Ref. jlj] in the unitary limit for an onsite 
potential. When an extended impurity potential is used 
with parameters similar to those obtained from the ex- 
act diagonalization calculations [|l8|, 2% impurities lead 
to the smearing of xoCQj ^) inducing spectral weight be- 
low wq- However, for 0.5% impurities a negligible effect 
is found. 

In the third section, the effects of the processes where 
the spin fluctuations scatter from the impurities with fi- 
nite momentum transfers were taken into account. It has 
been shown that the scatterings of the spin fluctuations 
with momentum transfers Q* 2ki? lead to a peak in 
Imx(Q,Li;) at cj ~ cjQ = 2|/i|. Here, the dependence of 
the impurity induced changes on the model parameters 
have been studied in order to make comparisons with the 
neutron scattering data. For instance, it was seen that 
the line shape of the peak depends on the model param- 



eters. When the parameters are such that the Stoner en- 
hancement is small, a dip is observed below the peak. On 
the other hand, when 1 — [/Rexo(Qi'^) is small, the dip 
is not observed and the 2k^ scatterings lead to a peak 
in Imx(Q,i^) by suppressing the damping of the RPA 
enhancement. In addition, in this case, the width of the 
peak increases with the broadening of the single-particle 
excitations. It has been also found that the impurity 
scatterings with finite momentum transfers could lead to 
spectral weight below wq. Hence, along with the results 
seen in the second section, the umklapp processes could 
also play a role in inducing the low frequency spectral 
weight observed in 2% Zn substituted YBa2Cu307 ||]. 

The general features of the impurity induced changes 
in Imx(Q,a;) calculated here appear to be in agree- 
ment with the experimental data on Zn substituted 
YBa2Cu307. However, in obtaining these results the 
Coulomb correlations were treated within RPA and 
Xo(Q, q, t^) has been calculated in the lowest order in the 
strength of the impurity potential. Furthermore, a static 
effective impurity potential was used. If the results pre- 
sented here are supported by higher order calculations, 
then it would mean that a contribution to the peak ob- 
served in Imx(Q, w) could arise from the magnetic chan- 
nel. Furthermore, in this case, the experimental data of 
Refs. |,§ would mean that a perturbation in the den- 
sity channel as Eq. (^5|) induces important changes in the 
antiferromagnetic response of YBa2Cu307. 
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